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1 Introduction 

In this work we find an asymptotic formula for a Toeplitz determinant with the symbol 
supported on an arc of the unit circle in the case when the symbol has Fisher-Hartwig 
singularities. 

Let f{6) be an integrable function on the unit circle. It was proved by Szego PP that if 
f{6) is positive and sufficiently smooth, namely, its derivative satisfies a Lipshitz condition, 
(these requirements have been later relaxed) then the Toeplitz determinant 

D4f) = det e-'^^-''^'f{e)de^ , j, k = 0,1 . . . ,n, 

has the asymptotic expansion 

n^oo, (1) 

where 



H{f) = exp (ln/)o, = exp fc(ln /)fc(ln 

k=l 

(In /), = ^ r e-'"' \nf{e)de, A; = 0, 1, . . . 
/vr Jo 



(2) 



This formula is not valid for f{6) with zeros or singularities. The asymptotics in the case 
when f[6) has zeros or root-type singularities were conjectured by Lenard i2j. Fisher and 
Hartwig and proved by Widom Namely, let 

20r 



fio) = m 11(2 - 2cos(^^ - e,)r = m n 



r=l 



r=l 



2 sin 



9-9r 



(3) 



where ip{6) satisfies the conditions of Szego's theorem and all > —1/2 (this is the most 
general case of real a^'s since it is a condition for existence of Fourier coefficients). Then 



2 sin 



(4) 
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where G{z) is Barnes' G- function. This was further generahzed for the case of f{6) with 
jumps by Basor Bottcher and Silbermann |H], and other authors. An incomplete hst 
includes the works [21^^01 • Note that the positivity condition for f{6) can be removed: the 
asymptotics (jH) are actually proved in |lj for a complex- valued f{9) with ip^O) ^ and the 
change of the argument over the closed circle A_.,r<9<n sugip^O) = 0. 

It is interesting how the asymptotics look like in case when f{6) is zero on an arc. 
This question was addressed by Widom in [TT]. Suppose that is supported on an arc 
a < (f < 27T — a, where it is positive, smooth enough, and symmetric: f{f) = /(27r — (f). 
Then it is proved in jTTj that 

~ r'H{Fr (nsin 1)"'^' ^^i^EiF), (5) 

with 7 = cos(q;/2) and 

Q 

F{6) = /(2 arccos(7 cos -)), 
the function /(v?) "extended" to the whole circle. Note that as is known 



where C'(^) is the derivative of Riemann's zeta function. Widom's proof was based on 
a theorem of Hirschman for a Toeplitz-like determinant, where in the matrix elements the 
Fourier coefficients of f{6) on the circle, i.e. the coefficients w.r.t. the system z^, are replaced 
by the coefficients of a function on the interval [—1, 1] w.r.t. the Legendre polynomials. 

Just like this formula breaks down if /(v?) has zeros or singularities on the arc. 
Resolution of this question is presented here. Let 

R 

/((/^) = V'lv^) n(2-2cos((/^-(/^,))"% a, >-l/2, = a, ^R = 2n-a, (6) 

r=l 

be supported on an arc a < (p < 27i — a, f{ip) = f{27i — ip) (in particular, ip^ = 2n — (pfj^i^r), 
and let ip{ip) satisfy the conditions of Szego's theorem. Then 



-^-.(/) ' ("-f)"'''^--^*--')^^ 



^ 6-^ + 1)2 / sin2(<^,/2) 



2 sm 



(7) 

where 

"^{9) = ^/'(2 arccos(7 cos -)). 
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In particular, when zeros and singularities are absent, we reproduce the asymptotics (0). 

For technical reasons, we need to replace ipif) by an analytic function. Smoothness of 
ip{(p) implies existence of a trigonometric polynomial ipoi^p) such that fi{{p) = ip{ip) —ipo^ip) = 
0{1/C) uniformly in on the arc for some constant C{n) which can be chosen arbitrary 
large. Let, e.g., C{n) = n". Then a simple analysis of the r.h.s. of ((Zj) (using smoothness 
and positivity of ip{(p)) and the l.h.s. of ((Tj) (using Hadamard's inequality) shows that it 
is sufficient to prove the theorem when ip is replaced by ipQ. To simplify notation, we just 
assume in what follows that ip{tf) is analytic. (Note that we have to bear in mind that ipif) 
now depends on n in a particular way, as it now stands for an analytic approximation to the 
original 

To obtain ((Tj), we use (jlj) to evaluate D„_i(F), where 

m = ,^^^^f^^y^. ^ ' = /(2arccos(7cos^)), (8) 

a -72cos2(0/2) 2 



is a symbol on the whole circle, and the following relation 



Dn-iif) _.„^+. Pn(7-^) .,n^ V/(^)(l + sm(a/2))"2^"T(l + 2«0 



Dn-i{F) Pn{l) A{F)^7msm{a/2) 

where 

= (27)^°V(^) 

n.=i,«(2-2cos(v.-^,))-' ^ ' 

sin(a/2) InF(g) 
^^^^^^-^^io l-ycos2(^/2)^^' 

and Pk{x) = + ■ ■ ■, k = 0,1, . . . are the monic polynomials orthogonal on [—1, 1] w.r.t. 
the weight function 



/(2arccos7x) 

= — n 2 2 • (12) 
V 1 — 7^a;^ 



Pn{x)Pmix)w{x)dx = K^mn, K > 0, m,n = 0,1, . . . (13) 

In our case of a symmetric weight {w{x) = w{—x)), these polynomials are simply related 
(see [121 Uni CH 113) to the polynomials ^^{^,(3) = + ■ ■ ■, k = 0,1, . . . orthogonal w.r.t. 
g^ip) = /(2 arccos(7 cos ^/ cos|)) on an arc: 

1 /•27r-/3 

— / ^niz)^rniz)g{(p)dip = hnSmn, ^ = c'^, /l„ > 0, TTl, U = , 1, . . . (14) 

ZTT J (3 

The first equation in © follows from the connection between 1) Toeplitz matrices, 2) 
orthogonal polynomials on an arc of the unit circle, and 3) those on a segment of the real 
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axis. To obtain the second equation, we need to analyze polynomials Pk{z) asymptotically. 
This is partly achieved with the help of an approach based on a matrix Riemann-Hilbert 
problem, a new method which turned out to be very effective for analysis of asymptotic 
problems fU^ IT7| HT^ . Most Riemann-Hilbert analysis we need is contained in fOl 1211 where 
the asymptotics were found for polynomials orthogonal on [—1, 1] with a weight (1 — a;)"(l + 
xYh{x), h{x) is positive and real analytic. We get the main term in the asymptotics of -Pn(l) 
by generalizing a small part of 120] to the weight p2|) . The main term for P„(7~^) can 
be obtained much easier, from the well-known Szego asymptotics (^21, Thm.12.1.2) which 
apply here because 7"^ lies in the exterior of the orthogonality interval [—1,1]. As will 
be clear from the proof, equations © also hold when f{9) has jumps in addition to zeros 
and root singularities (see 0). Therefore the result of Basor j3] can also be extended to 
symmetric symbols on an arc in the same way as Q. 



2 The ratio D^_i{f) / Dn-i{F). 

In this section we obtain equations ©. Due to the above mentioned connection between 
polynomials orthogonal on an interval of the real axis and on an arc of the unit circle one 
can write the following relations (Lemma 3.1. of jT5]): 

D^.,U) = 2"("-^)^P.(l/7) n hj, f{^) = w{r' cos f ) sin ^; (15) 

I)„_i(F) = 2"("-i)4^"(l) n^^' F(^) =^(cos^)sin^. (16) 

In the first of them the corresponding arc is a < y9 < 27r — a, and in the second one, the arc 
coincides with the whole circle < < 27r. Setting 7cos(^/2) = cos(<y9/2) and comparing 
the expressions for / and F we get (jS)). The ratio of D„_i's yields the first equation in Q. 
Note that we can also rewrite it as follows (using, e.g., formula (14) of [T3]): 

^n-i(/) ^ 2$„(l,a)7-'+" 
Dn-iiF) e-W2$^(ei",a) + C.C.' 

Thus, to obtain the second (asymptotic) equation in Q, we can either consider the asymp- 
totics of polynomials on an arc or on an interval. We choose the second option since in that 
case we can use the results of Kuijlaars, McLauphlin, Van Assche, and Vanlessen |l9j instead 
of reformulating the Riemann-Hilbert solution for polynomials on an arc. As indicated in the 
introduction, we need these methods to obtain the main term in the asymptotics of Pn(l)- 
In what follows, we mostly dwell on the differences of our case from |Tn| and refer the reader 
to that work for other details. 



Consider the 2x2 matrix 
Y{z) = ' 



27rm;_iP„_i(2;) /_i Pn-i{x)w{x 
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where Kn = h'^^"^ is the leading coefficient of the orthonormal polynomials. Y(z) is the 
unique solution of the following Riemann-Hilbert problem (the proof is almost the same as 
in ^H], the difference being additional singularities inside (—1, 1)): 



(a) Y{z) is analytic for z G C \ [—1, 1]. 

(b) For X G (—1, 1) \ UrXr, Xr = ■j^^ cos(v5r/2), Y has continuous boundary values Y^{x) 
as z approaches x from above, and Y^{x), from below. They are related by the jump 
condition 

r+(x)=F_(x)Q ""f^y Xei-lA)\UrXr. (19) 

(c) Y{z) has the following asymptotic behavior at infinity: 

F(z)= (j + oQ)) ("J asz-.oo. (20) 



(d) Near the points 1 = xi, —1 = xr, Xr 



Y(z) = 



O 
O 



2ar ; ' 



1 1^ — XrP"'' 

1 log \ z — Xr 
1 \og\z — Xr 



o 



1 1 
1 1 



if Ur < 0, 
if ttr = 0, 
if ttr > 0, 



(21) 



as Z ^ Xr, z & C \ [—1, 1]. 



Here the points Xr correspond to ipr under the mapping x = 7 ^cos{ip/2). Moreover, 
recalling (0), (fT2j) . we have 



w(x) 



ip{2 arccos7x) 



R 



X 



ICLr 



(22) 



The proof of (j2H) is similar to that in jl9j . 

Now we find the asymptotics of the matrix element Y\\{^z^ in the neighbourhood of z = 1. 
To do this, we apply a series of transformations to the Riemann-Hilbert problem for Y{^z\ 
as is usual in this method. These transformations Y ^ T ^ S are the same as in but 
the contour where the jump for S'(z) occurs is different. It consists of i? — 1 lenses which 
touch at zeros and singularities Xr of the weight Wq[x) (see the figure). In ^U] there was only 
one lens with the endpoints —1 and 1. We now consider the neighbourhoods of the points 
Xr and the region outside of them separately. In the outside region, a parametrix Ni^z^ is 
given by the same expression as in JH] where we substitute our weight ^0(2;). It coincides 
with the main asymptotic term for ^(z) obtained from Szego's Theorem 12.1.2 in jl2j. In 
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Figure 1: The contour for the Riemann-Hilbert problem is sketched in the case of i? = 3. 



the neighbourhood of each point x^, one could solve the local Riemann-Hilbert problem such 
that it fits N{z) on the neighbourhood's boundary. Of those solutions, we need only the 
one in the neighbourhood of xi = 1. Indeed, it is easy to show again like in Theorems 7.8, 
7.9 of 121] that corrections from neighbourhoods of the other points Xr contribute at most 
to the next-leading term in the asymptotics near x = 1. Therefore the leading term near 
X = 1 coincides exactly with that for the problem considered in l2Uj. It is given by the 
formula (2.23) and the one after (3.7) of [20]. Taking there the limit x — > 1, we obtain the 
asymptotics of Yii{l) and finally, recalling (fTHjl . get the relation: 



where 



"^^ ~ 2" V 1 + sin(a/2) \/ v{l) 2^-^T{l + 2a,) ' ^ ' 

, „ 1 /"^ In f (2 arccos7x) , 1 /"^'r 

and f (x) = w(x)/(l — x^)^"^ 

One can verify that there is no problem in the fact that ijj{(p) (and hence w{x)) depends 
on n (in a special way as it approximates some smooth function). 

We now turn to calculation of the asymptotics for P„(7~^). Here it suffices to use an 
old result of Szego. The leading term for asymptotics of ori/ionorma/ polynomials Pn{x) for 
X outside [—1,1] and a wide class of weights on [—1, 1] (which includes our weight (|221)) is 
given by Theorem 12.1.2 in ^21- Namely, 



1 

~ /7^^/ — ?T' a; = -(z + z-^), (25) 
v27rPU^ 2 



\nV{z) = — ln(w(cost)| sint|)- —dt. (26) 

47r J-TT 1 — ze ** 

For X = 1/7 = l/cos(a/2) we have z = {1 + sin(a/2))/7 = zq. Substituting ^ into P^ . 
separating the integral into 2 parts, one over [0, tt], the other over [— vr, 0], and changing the 
variables 6 \—>- —6 in the second one, we get 

, ^/ ix ^ r f sin 61 \ 2sin(a;/2) ,^ , ^/r^x 

lnV{z^') = — In F{29) == ^-^d9 = In A F + J, 

47r JO \ V 1 — 7 cos"^ 6* y 1 — 7 cos 

where 

_ sin(a/2) lnF(2g) _ sin(a/2) r^^ InFjO) 
^ ' 2-K Jo 1-7COS0 47r io 1 - 7^ cos2(^/2) 
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(we used the symmetry F{0) = F{2n — 6) to get the last equation) and 

sin(a/2) /•2- Isin^l 69 



4:71 Jo \/l — 7^ cos^ 9 1 — 7 cos 9 
The last integral can be evaluated explicitly. We have for |a;|>l,0<a<l 

Vx2^ /•2-ln(l-acos2^)^^ , s + ^i^^ 

' -a9 = In ■ (28) 



An Jo X - cos 9 - 1 + xv^l - a 

This integral is calculated by differentiating w.r.t. a and then computing residues in the 
complex plane of e*^. Taking here a — > 1 and setting x = I/7 gives part of the integral (j27j) 
with In I sin 6*1. The part with ln(l — 7^ cos^ 6') is obtained by setting a = 7^, x = I/7. Thus 
we get J = — In a/2, and therefore 

\nV{z,')=\n^. (29) 

The main term in the asymptotics of the leading coefficient k„ of Pn{x) = KnX"' + ■ ■ ■ is given 
by Theorem 12.7.1 of A calculation like the one just carried out yields 



2" l + sin(a/2) 



where is given by (j^ . Putting (P^. ^2^ . and (jSHI) together, we finally obtain 



K„ V 27 / V 1 + sm(a/2) A(F) 
Equations and (p?T|l give the asymptotics 0- 



3 Asymptotics for Toeplitz determinants 

We are now ready to obtain the asymptotics (|7j). An important role will be played by the 
map between the arc and the circle given by cos(y9/2) = 7cos(6'/2). Using the symmetry of 
the symbol f{ip) = /(27r — ip), and hence yj^ = 27r — ipji^i_r, we can write (jH)) as {9rg = vr) 



[R/2] 



F{9) = /(^) = ^(y.) n 



r=l 



4 sm sm 



2q,. 



2 COS 



2ar 



[R/2] 

^{^) n (4 

r=l 



2 2 V^r 

COS COS — 

2 2 



2ar 



2 cos 



2ar 



(32) 



R 



2 sin 



2ar 



r=l 
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Consider the function F (jH)). Vanishing of /(v?) on the arc —a < ip < a is "replaced" in 
9 variable by a singularity of F at ^ = given by the factor 2sin(^/2). If f{<p) also has 
singularities at the endpoints of the arc, the corresponding factor becomes (2 sin(0/2))^^^"^ 
We see that Widom's formula (jH) holds for F . The role of if) is now played by 
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2^1 -72cos2(^/2) 

tti in the formula is replaced by 1/2 + 2ai, i?, by i? — 1 (as the ends of the arc merge in 9 
variable). Let us now simplify this expression for Dn-i{F). We have 



In 2W 1 — 7^ cos^ 



These integrals are calculated as 



9 



1 / l-sin(a/2) \ 
'm \ 1 ) 



2\k\ 



k = 0. 



In 2 



ln(l + sin(a/2)), 
. In particular, if we set a = for k = 0, we get 

9 — 9r 



(33) 



sm 



0. 



The first obvious application of these relations is the expression 



H{F) = H 



7 



2Er"'-^(6') 



7 



2Er°'-i^(^) 



2^1 -72 cos2(0/2) J 1 + sin(a/2) 



(34) 



According to (g)), we now have to calculate E(^(6')(l - 72 cos2(6'/2))~^/2)_ (^^ote that 
E{Cv{9)) = E{v{9)), where C is any nonzero constant.) Using (j33|l and performing sum- 
mation of geometric series, we obtain: 



EMM*W(l-7'cos^?)-i/n 



fc=i 
1 

2"7r Jo 



In^(^) ^ ^ 



r 



sinf)2e-^^"" 4 



k=l 

rf^ - - In 1 1 - 7~^ 



(35) 



a 



sm 



The imaginary part of the integral here vanishes while the real part can be written as 



--(lnvl>(^))o + lnA(vl/), 

where A(-) is defined in (fTT|) . 

Exponentiating (jHSj) and multipliyng the result with (2sin(Q;/2)/(l + sin(a/2))^/2^ 
finish calculation of E: 



2sin(a/2) 
^ 1 + sin(a/2) 



E 



^>{9) 



a 



72 cos2(0/2)^ 



i?(v^)4i£LsinV^ 



(36) 



Separating the contribution of the point 6'i = 0, we can now rewrite (0} as follows: 



Dn-l{F) 



l + sin(a/2) / 



l^(H{m) 2sin((^,/2) G{ar + l 



,2 -R-1 

1 + sin(a/2) ; i-^, G{2ar + 1) 



n 



2 sm — 
2 



^^(0) 



2 G(2 + 4ai 



n 



2 sin 



— 0^ 



(37) 



From (UnD, (122) we obtain 

/»=2^"^f ^^"'-^w n 

The ratio Q for this /(v^) can be written as 

Dn-l{f) 



2 sin 



6 — Or 



2Qr 



Dn-l{F) 

where 
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i:.-'-V^nf=2'(2sin(0,/2))"- 



a\"24"ir(l + 2ail 



1 + sin ■ 

A(^)iry'7msin(«/2) 

ln|27sin((0 -0,.)/2)| 



n 



2ai 



47r 



1 -72cos2(0/2) 



d9. 



(38) 



(39) 



Using the symmetry of ^^'s we regroup the addends in the above expression and write 

6 — 6r . 6 -\- 6r 



In 



4 sin 



sm ■ 



In |1 - 26cos6l| - ln6, l/6 = 2cos6'^ 



2 2 

Now the integrals can be computed as before, and we get 
Putting this together with (jHHj) and (jHTj) . we obtain 



(40) 



2^"ir(l + 2ai)n 



v/^G(2 + 4ai)' 
^ /^i/(^)V-^j^^ gK + l)^^ (2sin(g./2))-^°^°- 



7/;(y.,)y VG'(2ar + l) 



2 sin 



^^r — ds 



(41) 



Finally, a simple but cumbersome calculation which uses the symmetry of O^s (cf. 1)32 
gives 



n 



2 sin 



ER-1 
2 "-"^-Z^ 



a2 

2 



n 



2 sm 



Hi 



sin^(v9r/2) 



(42) 



r = 2 



— 7 2 C0s2((y9r/2)_ 



Substituting this into the previous formula completes the derivation of ((Tj). 
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